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recently constructed by two of the authors as an outcome of a relation among FZZT 
branes. In the double scaling limit, the IMM is described by the (2, 2p+ 1) minimal 
Liouville gravity. These matrix operators are shown to create a boundary with 
matter boundary conditions given by the Cardy states. We also demonstrate a 
recursion relation among the matrix disc correlator with two different boundaries. 
This construction is then extended to the two matrix model and the disc correlator 
with two boundaries is compared with the Liouville boundary two point functions. 
In addition, the realization within the matrix model of several symmetries among 
FZZT branes is discussed. 
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1. Introduction 

Since their introduction in the middle 80 's in the context of string theory |1[, 0, ^ and 
lattice statistical models, 0, ^ § , the role played by the boundary conditions (BC) in 
2D conformal field theories (CFT) has been widely uncovered 0, [^. However, several 
intriguing questions remain, the most challenging involving the interplay between 
boundaries and renormalization group (RG) flows. This includes the description 
of RG flows driven by perturbations of conformal boundary conditions 0, and the 
matching of boundary conditions between two CFTs related under a bulk RG flow 
1^. Both attempts to tackle these problems, which rely respectively on integrability 



and perturbative RG analysis, make use of the ^f-function to identify the boundary 
conditions. Unfortunately, the results obtained are restricted to the thermal flows 
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of the unitary {p,p + 1) minimal models, more precisely to the tricritical to critical 
Ising flow in 0, and to a pertubative expansion for large p in Very recently, 



these results have been extended to any p under plausible assumptions, using the 



staircase model 11 



We focus here on a very promising alternative approach based on matrix models. 
Matrix models can be used to construct statistical models defined on a dynamical 
lattice and allow to take an explicit continuum limit. In this double scaling limit, 
the critical points are described by a CFT, referred as "matter", coupled to a two 
dimensional quantum metric. In the conformal gauge, the metric field reduces to 
a scalar field that obeys the Liouville action, and the theory is called "Liouville 
gravity" . The matter CFT is nothing but the theory that describes the scaling limit 
of the statistical model on a fixed lattice at its critical point. The study of the 
Liouville gravity correlation functions allows to recover several features of the bare 



matter CFT by making use of the KPZ relations ||T2|, |13| |14| between the central 
charges and operator dimensions of matter and Liouville theories. Since the matrix 
model description is also valid away from the critical point, it reveals a profitable 
way to understand both bulk and boundary RG flows. This approach has been 
successfully applied to the lattice 0{n) model whose continuum limit consists in 
a set of non-rational CFT which interpolates between the unitary minimal series. 
There, the matrix model gave a description of the thermal bulk [[l^ and boundary 



Iq] flows, as well as the boundary anisotropic flow |jT6|, |T^, |l^. The goal of the 
present article is to pave the way for a similar construction within the two matrix 
model (2MM). Taking a suitable double scaling limit, this model is known to describe 
a (p, q) minimal model coupled to the Liouville field, sometimes referred as (p, q) 
minimal Liouville gravity. 

In the Liouville gravity framework, the boundary condition changing operators 
(or bcc operators) ipa of the matter theory are dressed by the Liouville field ip: 
the operators are multiplied by a boundary vertex operator of the Liouville theory, 
Ba = V'ae^"'^. In order to get a diffeomorphism invariant operator, the coordinates 
should also be integrated over the boundary of the worldsheet. In the following, we 
will consider the matrix model in the planar limit, thus restricting to worldsheets 
with the topology of a disc. The Liouville charge (3a is related to the dimension 
ha of the matter boundary operator through the KPZ formula. In the case of the 
{p,q) minimal models {p < q for definiteness) , the bcc operators fill the Kac table: 
a = {k, I) with 1 < k < p and 1 < / < g. In addition, ipk,! and ipp-k,q-i should be 
identified since they wear the same scaling dimensions 

hk,i = hp_k,,-i = \ [{k/h - Ihf - (1/6 - hf] , I? = p/q . (1.1) 
The KPZ equation is quadratic in (3a, leading to two different solutions, or "dress- 
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mgs ' 



Both dressings can be realized in the 2MM and we will see in particular that the 
symmetry (A;, /) — {p — k,q — l) involve an exchange of the dressings for the boundary 
two point functions. However, when not mentioned explicitly, the dressing will be 
supposed to satisfy the Seiberg's bound ||19|, (3 = f3~^. We will come back to this 
subtle point below. 

The boundary states must specify the boundary conditions for the matter and 
the Liouville fields. Called FZZT branes they are constructed as a tensor 

product of a Cardy brane \k,l > with a Liouville boundary state \s >. The boundary 
parameter s is related to the cosmological constant ^(s) which control the length of 
the boundary. The presence of conformal Killing vectors allow to fix the position 
of three operators on the boundary. As a consequence of the absence of interaction 
terms between the Liouville and matter fields in the action, both contributions to 
the boundary one, two and three point functions factorize. The Liouville correlator 
is the only part that contains a dependency in the cosmological constants. Most of 
the results of this paper exploits this property of the boundary two point functions. 
More precisely, the Liouville momentum / can be identified by comparing the 
cosmological constant dependence of the matrix model correlators in the continuum 
limit and the expression for the Liouville boundary two point function which has 



been derived in ||20|. This identification provide a way to determine the dimension 
of the bcc operator inserted between two given boundaries. 

The first step of this ambitious program to study boundary RG fiows is to con- 
struct the equivalent of the Liouville gravity bcc operators within the matrix model. 
This has been recently achieved in the one matrix model (IMM) p2| and will be 
reviewed in the first section. This matrix model can be tuned to give a description 
of (2, 2p + 1) minimal Liouville gravity in the double scaling limit. The aim of the 
present article is to discuss further this construction and to extend it to the 2MM 
case (section 2) in order to recover the bcc operators of general {p, q) minimal Li- 
ouville gravity. The strategy we employed is first to introduce the matrix operator 
that creates a boundary inside a single boundary correlator, motivated by a decom- 
position relation among FZZT branes Then, we turn to two boundary matrix 
correlator and determine the bcc operators inserted from the identification with the 
boundary Liouville two point function, as explained above. We eventually deduce 
the matter boundary conditions from a CFT argument. For a matter of clarity, 
most of the calculations are confined to the appendix. It contains in particular the 
demonstration of the two main results: the recursion relation ( 2.10 ) for the IMM disc 
correlators with two general boundaries already mentioned in p2[ but not proven, 
and the identification of the 2MM correlation functions ( p.l3| ) with the correspond- 
ing Liouville boundary two point function. Additional remarks and perspectives are 
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gathered in the conclusion. 



2. One matrix model 

The IMM partition function is defined as an integral over the hermitian matrix M 
of size N X N, with an action given by a polynomial potential V{My, 

'-ti V{M) 



D[M] e-^''""^''''. (2.1) 

The free energy A^~^ log (Z) is a generating function of discretized closed surfaces 
where the weight of the polygons is given by the polynomial coefficients of the corre- 
sponding degree. In the large expansion of log (Z), each surface carries an overall 
factor K~^N^ depending on its topology through the Euler characteristic x? ^ind 
on its area A with a cosmological constant = 1/g. The planar limit consists in 
sending toward infinity, keeping only the surfaces with the sphere topology. The 
main quantity of the model is the spectral density associated to the matrix M, it is 
usually given through its Stieltjes transform, the resolvent W{x): 

W{x) = ^/ti = ^^^^ + u{x) . (2.2) 



\ x-M, 

In the planar limit, the resolvent describes a sum over discretized surfaces with the 
disc topology and one point marked on the boundary. The boundary cosmological 
constant x controls the length of the boundary, as can be seen from the a; — t- oo 
expansion. When the number of eigenvalues goes toward infinity, the support of 
the density condensate into a continuous subset of the real line and the resolvent 
develops a branch cut. For our purpose, we concentrate on the single cut case where 
the support is an interval (for a review, see |2^). 

The continuum limit is achieved by sending the bulk and boundary cosmological 
constants to their critical values k* and x*. Note that by shifting the matrix, the 
potential can always be chosen such that the critical value of x is zero. The minimal 
(2, 2p+l) Liouville gravity is obtained after tuning the coefficients of the potential up 
to the degree p in the double scaling limit. At the critical point, the mean area and 
boundary length diverge and we define the appropriate renormalized cosmological 
constants e^yU = k — n* and = x, with 6^ = 2/{2p + 1). The singular part of the 
resolvent^ e^^^'^uj{^) = u{x) defined in ( |2.2| ) is identified to the one point function of 



the boundary operator Bi i on the disc. The renormalized quantity ^ corresponds to 
the Liouville boundary cosmological constant, sometimes denoted /ib in the literature 



^In the following we consider (2.1) as a formal matrix model, i.e. the matrix integrals are seen 
as formal series which terms are computed using the gaussian kernel. 

^The non critical (or non universal) part, which is always polynomial, can be related to degen- 
erate configurations and shall be discarded in the continuum limit, see [ESl. 
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pO[] . In order to resolve the branch cut of the resolvent on ] — cxd, —u], ^ is usually 



parameterized by the boundary parameter s, 

^{s) = u cosh {Tihs) , = u^^^^ cosh (775/6) . (2.3) 

where u = a/ /u/ sm{iTb'^). We will assume that the matter carries (1, 1) BC on the 
disc, which will be confirmed by the consistency of the construction below. 

2.1 Construction of the matrix boundaries 



In 1^ , non-trivial matrix boundaries were introduced using an additional vector field 
which couple to the matrix M. We use here another approach based on a relation 
among FZZT branes first noticed in |]2B|, and derived in from the boundary 
ground ring technique. This relation decomposes the FZZT brane |s; 1, Z > as a sum 
over FZZT branes 1, 1 > with trivial matter BC and shifted boundary parameter 
Sn = s + inb, 

l-l 

\s;l,l>= Yl k";l,l> (2.4) 

n=-(«-l),2 

where the index n runs from — (/ — 1) to (/ — 1) with an increment of two. The 
relation ( p.4| ) trivially extends to the disc partition functions of the Liouville gravity 
Zlg{s] 1,1). The resolvent is the derivative of the disc partition function Zi^q{s; 1, 1) 
with respect to the cosmological constant. Consequently, the partition function 
Zlg(s;1,1) is obtained by integration of (|2.2|), i.e. by taking the double scaling 
limit of ;^ (tr log(x — M)). Applying the linear relation (|2.4|) to the logarithmic 
correlator leads to regard the following quantity as a candidate for the more general 
Zlg(s; 1, /) disc partition functions, 

l-l 

Zi{{x}i) = ^ ^ (tr log (a;„ - M)) , = cosh (vrfosn) . (2.5) 

n=-{l-l),2 



In this way, we recover the partition function considered in where was defined 
the matrix operator Fi{M) as 

l-l 

Zi{{x}i) = ^ (tr logF^M)) , Fi{M) = J] (x„ - M) . (2.6) 

n=-(«-l),2 

We will show that the operator Fi creates a boundary with cosmological constant 
X = eu cosh [Trbs) and matter boundary conditions (1,0- First, we notice that the 
derivative with respect to x of the disc partition function Zi{{x}i) is proportional to 
the resolvent, and thus corresponds to the Bi i boundary one point function, 

Oi{{xh) = j-d. (tr logFKM)) = {-iy-'[l]Mx) + (n.u.) , [/]+ = . 

iV sm (tto^) 

(2.7) 
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Here (n.u.) refers to the non universal polynomial part of the correlator. This con- 
firms the identification of the cosmological constant x. The assertion concerning the 
(1,1) matter boundary conditions is supported by taking the ratio^ of the singular 
part of Oi and of the resolvent |22||. We thus obtain the factor (— which is 
exactly the ratio of the matter parts / ('^i,i)(ii)- Hence, the correlator Oi 

corresponds in the double scaling limit to the Bi i boundary one point function with 
^(s) cosmological constant and (1,/) matter BC. This elucidates the role played by 
Fi at the level of the disc with a single boundary. 

In obtaining ( p.7|) , it is crucial that the resolvent take the same value over the 
shifted cosmological constants Xn, = {—iy~^u{x). This condition restricts to 

at most p different solutions for x„ so that / < p. Note that, as supposed by the 
identification of the boundary conditions (1, /) and (1, 2p + 1 — Z), the result ( p.7[ ) is 
invariant under the exchange I ^ 2p + 1 — I. 



2.2 Matrix correlators and Liouville boundary two point functions 

We investigate here the role played by the matrix operator Fi at the level of a disc 
with two boundaries. More precisely, we present a method to determine the matter 
boundary conditions associated to the boundary created by Fi, based on the study 
of the Liouville boundary two point function. This analysis will be extended to 
the 2MM case in the next section. The idea is to consider the following correlator, 
interpreted as a disc with two boundaries constructed with respectively Fi and Fi, 

In the continuum limit, the singular part of Oi is proportional to the Liouville bound- 



ary two point function of Bi i operators, as proven in |]22[. The matrix operator Fi 
creates a boundary with BC (1, 1) for the matter field. The matter bcc operator must 
belong to the fusion algebra of the Verma module associated with the two nearby 
boundary conditions. Any BC fused with the identity reproduces itself, thus, on the 
boundary created by Fi, the matter field BC wear the same label as the boundary 
operator: (1, 1). This identification is compatible with the symmetry I 2p + 1 — I, 
provided we exchange the dressings by the Liouville field. We will come back to this 
point in the next section, dealing with a slightly more general case within the 2MM. 

Once the role of the Fi matrix operator has been elucidated, it is possible to 
understand the general two boundary matrix correlators 



■'The ratio allows to neglect the contributions of leg factors and coordinates. Moreover, as the 
Liouville part does not depend on the matter BC but only on the dimension of the matter boundary 
operators, it also simplifies from this ratio. 
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The correlator Oi^k corresponds to a disc with (1,/) and {l,k) matter BC in the 
continuum hmit. These matter BC allows several bcc operator insertion, namely 
ipi^n where n = 1 + \l — k\ ■ ■ ■ rimax, 2 and rimax = niin{l + l + k,4:p+l — I — k}. 
In order to determine the bcc operators involved in 0;,^, we employ the following 
recursion relation, 

[k]+Oi,k+ii{x}i, {y}k+i) = [/]+Om,fc({x},+i, {y}k) + (n.u.) . (2.10) 



This relation was already mentioned in but not proven. Its demonstration 



is one of the main results of this paper, it can be found in the appendix (|A.3|) . 
Using the relation ( p.lO|) to decrease the index / up to one, we eventually come 



back to the known correlator Oi-^-k-i defined in ( p.8|) . The identification with the 
Liouville boundary two point function leads to determine the bcc operator 
We notice that this operator depends only on the sum / + k, implying that many 
different boundary conditions lead to the same bcc operator. This is the meaning 
of the recursion ( p. 101 ): both Liouville part of the gravitational correlators are the 



same whereas the factor [/]+/[/c]+ reflects the different matter parts. Unfortunately, 
it is not known yet how to construct the matrix model correlator corresponding to 
the insertion of the other bcc operators between (1, /) and (1, k) matter BC. This is 
one of the major open questions we hope to answer in a near future. 

Up to now, we have demonstrated that the matrix operator Fi can be used to 
replace Fi in the case of a disc with one or two boundaries, thus promoting (1, 1) 
to (1,0 matter BC. We believe that this can be generalized to any correlator of any 
genus, with an arbitrary number of boundaries, possibly disconnected. Moreover, the 
relation ( |A.1| ) employed in the proofs above allow to decompose these Fi boundaries 
as a sum over Fi boundaries with shifted cosmological constants, similarly to ( |2.4| ). 
All these correlators involving only Fi can be evaluated recursively using the method 



developed in [^. The main complexity in the computation of the general correlator 



lies in the evaluation of the sum over shifted boundary parameters. 



3. Two matrix model 



The two matrix model was first introduced in p9| in order to describe the 2D Ising 
model on a fluctuating lattice. It was subsequently generalized to arbitrary poly- 
nomial potential. The partition function is defined as an integral over two N x N 
hermitian matrices X and Y with potentials V{X) and U{Y), 



dXdY exp 



tr (V(X) + U(Y) ~ XY 

9 ^ 



(3.1) 



Each matrix is associated to a site wearing a different color and the XY propagator 
is orthogonal to an edge of the interface between the two colors. The potentials 
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generate propagators and vertices relating sites of the same color, similarly to the 
IMM. In the case of a quadratic potential U{Y), the matrix Y can be integrated 
out and we recover the IMM. The main quantities are the two eigenvalue densities 
associated to the matrices X and Y, so that we need to introduce the resolvents'^, 

These matrix correlators describe a disc with a boundary of a single color, X or Y, 
in the presence of one marked point. ^ Again, we restrict ourselves to the planar limit 
and the resolvents are multivalued functions with a branch cut on the support of 
the densities. We denote their respective singular part ijj{x) = V'{x) + W{x) and 
uj{y) = U'iy) + W{y)^ and the bare (bulk) cosmological constant is = 1/g. 

In the double scaling limit, the coefficients of the potentials V and U are tuned 
up to degree q and p respectively, in order to achieve the (p, q) critical behavior 
{p < q for definiteness) |^. Using a shift of the matrices, we set the critical values 
of the boundary cosmological constant to be zero. We then define the renormalized 
cosmological constants e^^/x = k, — k,*, e^^ = x and y = e'^(. The dependence 
in the cut-off as been slightly modified from the IMM in order to emphasize the 
symmetry between p and q that exchanges the two matrices X and Y. The boundary 
cosmological constants are conveniently parameterized as^ 

^(s) = 2uP cosh (ns/q) , C{t) = cosh (nt/p) , = p/q . (3.3) 

The resolvent u{x) is identified in the continuum limit to the Bi i boundary one point 
function of the minimal Liouville gravity, with ^(s) boundary cosmological constant 
and (1, 1) matter BC, in agreement with the results of the first section. The Y matrix 
resolvent also corresponds to the minimal Liouville gravity Bi^i one point function, 
but the boundary wears the dual cosmological constant C('S), sometimes denoted /is 
in the literature. The identification of the matter BC corresponding to u}{y) will be 
made below. Both resolvents can be easily expressed using the boundary parameters 
s and t, 

w(0 = 2m« cosh (tis/p) , u{C) = 2vF cosh (vrt/g) . (3.4) 

Most of the information concerning the 2MM can be encoded into a single object, 
the spectral curve defined by 

E{x, y) = {x- U\y)){y\x) - y) - P(x, y)^g, (3.5) 

^For later convenience we included a minus sign in the definition of the resolvents. 

^See | |30| ] for the amplitudes describing cylinder and higher disconnected loop-boundaries with 
ordinary boundary conditions. 

^Again, we slightly modify here the previous parameterization in order to emphasize the sym- 
metry. The boundary parameters of IMM and 2MM are thus related by 6simm — S2mm/<1, 
simm/& = S2mm/p- The variable C(i) is identifed with a dual cosmological constant, thus the 
role played by 6 and 1/6 are exchanged, ^imm/^ = i2MM/p, btiuu — ^2Mm/9- 
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with the polynomial 

^'^^A^\ x-X y-Y I ^ ' 

Since E{x,y) is a polynomial, the equation E{x,y) = defines an algebraic curve 
in C X C. As explained in the appendix ( p.l|) , this curve is equivalent to the con- 
dition X = u){y), or by symmetry y = 00 (x), that characterizes the distribution of 
eigenvalues. In the continuum limit, using the parameterization ( p.3| ), its expression 
is remarkably simple (see appendix (|B.1D for an explicit derivation). 



E{x, y) = eP%P« [cosh (tts) - cosh (vrt)] . (3.7) 

The spectral curve allows to determine recursively all the correlators of the model 
[0^ , ^ . More specifically, it will be related here to the expression of the mixed traces 
studied below. 

3.1 Study of the mixed trace 

The boundary operators of the 2MM were first considered in connection with the 
FZZT branes of Liouville gravity in [Q. In this paper, the author focus on the 



mixed trace correlation function, given by 

The derivation of this expression is reviewed in appendix (p.2|) . In the combinatorial 
interpretation of the 2MM, the mixed trace correlator corresponds to a disc with two 
boundaries of different color, each boundary is constructed using a different matrix 
(see e.g. [0). The singular part of M. is defined in the double scaling limit as 
eP'i^p^iM.{^X) = M.{x,y) — 1 which writes, using the parameterization (|3.3|), 

j^i^i C(0) u'P'^~^~'^ cosh (tts) — cosh (vrt) 

4 [cosh (vrs/g) — cosh (vrt/g)] [cosh (vrs/p) — cosh (vrt/p)] 

(3.9) 

In [^ , this expression was shown to coincide with the Liouville part of the gravita- 
tional -Bp-1,1 boundary two point function with boundary parameters s and t — ipq. 
From the supposition that the boundary created by X matrices has matter BC (1, 1), 
we deduce the matter BC {p — 1, 1) of the Y boundary. This claim will be checked 



at the level of the one point functions in the section |3]^ below. 
3.2 Analysis of the 2MM boundaries 

Since the 2MM can be seen as a generalization of the IMM, the later corresponding 
to the case p = 2, the analysis of the first section still remains valid. In particu- 
lar, the matrix operator Fi{X) leads to a boundary with (1,/) matter BC and ^{s) 
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cosmo logical constant. Note that these matter boundary conditions form a closed 
subset under fusion, 



(1,0 X (!,/')= ^max = min(l + / + r,2g-/-r-l). (3.10) 

n=l-\l-l'\ 



The decomposition of FZZT-branes (|2.4] ) found in ||26| can be generalized to the (p, q) 
minimal Liouville gravity, 

i-i fc-i 

\s]k,l>= ^ ^ |sm„;l,l>, Smn = s + imp + inq . (3.11) 

m=-{l-l),2 n=-(fc-l),2 

As in the one matrix model, the operator Fi{X) provides the case k = 1. The dual 
case / = 1 is introduced as a new matrix operator, 

fc-i 

My)= n (Vn-y), l/n(tn) = 2e^M^cosh(7rt„/p), tn = t + luq . (3.12) 

n=-{k-l),2 

In order to describe the boundary created by FkiY) in the double scaling limit, we 
will consider the following matrix correlators, 

M.(W,.{!/W = (tr^^). (3.13) 

The expression of M.i^k is worked out in the continuum limit in the appendix ( [B.3|) . 
This correlator is seen to reproduce the Liouville boundary two point function of 
Bp_k i operators with boundary parameters s and t — ipq, which is one of the main 
results of this paper. 

We first consider the case / = 1 where the left boundary has matter BC (1,1). 
The matter bcc operator inserted is ipp-k,i so that the right boundary, created by 
Fk{Y), must have {p — k,l) matter BC. This result is consistent with the study 
of ratios of Bi i one point functions (see section |3.3|) . We notice that this set of 
boundary conditions is also closed under fusion. 

We are now able to interpret the double scaling limit of the correlators Mi^k with 
any / as a disc with two boundaries with boundary parameters s and t — ipq related 
respectively to the cosmological constants ^ (usual) and ( (dual). The matter wears 
boundary conditions (1, /) on the left side and {p — k, 1) on the right side. Between 
both, the bcc operator Bp-k,i is inserted. This operator is the only one that can be 
inserted between (1, /) and {p — k, 1) matter BC. Varying k and /, we can realize all 
the operator dimensions belonging to the Kac table. 

It is also possible to consider the matrix correlators made of Fi{Y) and Fk{Y), 

o«(W<,{yW = (tr^^) (3.14) 
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The description of this quantity can be derived from the symmetry that exchanges 
p and q, reversing X and Y matrices. The correlator Ok,i thus represents a disc 
with {p — 1, 1) and [p — k, 1) matter BC and dual boundary cosmological constants. 
Between the two boundaries, the operator -B2p-fe-«-i,i is inserted. 



A priori, the relation ( |3.11| ) allows to construct boundaries with any matter BC 
{k,l). Unfortunately, a naive application of our method does not work, due to the 
independence of the cosmological constants ^{smn) in n (or similarly of C{smn) in 
m). This is not a problem in Liouville gravity because the Liouville BC is really 
characterized by the boundary parameter Smn, instead of the cosmological constant 
^(smn), due to the quantum corrections ||2T[. On the contrary, the 2MM correlators 
with (1, 1) matter BC depends on Smn only through ^(s^n)- 

To be more precise, the problem is well defined at the level of the spectral curve. 
There, the equation E{^,u) = gives q solutions in ^ at fixed u, and p solutions in 
u at fixed ^ : 



UJ„ 



2 cosh 
2 cosh 



— (s + imp) 
1 



— (s + inq) 
p 



m 



n 



-(g-l)---(g-l),2 
-(p-l)---(p-l),2. 



(3.15) 



We thus have determined q values of ^ and p values of uj such that E{^rn, ^n) = 0. It 
is then possible to choose / contiguous solutions for ^, i.e. living on adjacent sheets, 
and k solutions for co, centered on ^o, ^o, in order to construct a boundary with {k, I) 
matter BC. For instance, the boundary one point function will be given by 



/-I 



fc-i 



UJr. 



din 



\l+k 



m=-(Z-l),2 n=-(A:-l),2 



di 



sin (nb^'^k) 



sin (7r6='=^) 

(3.16) 

with ^ = 2m*' cosh (vrs/g) and u;(^) = 2m'' cosh(7rs/p). We verify that the quotient 
d^Zk^ijuj^i) corresponds to the matter boundary one point function contributions, 
expressed as a ratio of the transformation matrix of characters |3^ 



-1 



xl+k 



[MkV 



(3.17) 



In a similar construction, the disc with two boundaries, cosmological constants ^ 
and and matter BC respectively (/c, I) and (1, 1) is obtained from the following 
expression. 



V 



(e,c) 



Y\m=-(l-l'].2 [C 



-1 



(3.18) 



lm=-(Z-l),2 LS ^ V^-^; in 

It is easy to check that it indeed gives the Liouville boundary two point functions of 
Bk I operators. However, in the 2MM settings, only the "physical solution" ujq among 
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the p possibilities for the resolvent is selected, due to the asymptotic requirement 
before the continuum limit, W{x) ~ Nevertheless, Zki and T>{^k,i),(i,i) might 

be obtained as a proper matrix model correlator. 

3.3 Symmetries among FZZT branes 

In this section, we investigate the realization of various symmetries of the FZZT 
branes within the matrix model. The following discussion can be applied to both 
one and two matrix models. 

Let us first concentrate on the equivalence between matter boundary conditions 
(fc, /) and {p — k,q — I). In the appendix, we show that the correlation function Aik,i 
can be identified either to the Liouville two point function with a charge Pp-k,i, or to 
the inverse of the Liouville two point function with a charge Pk,q-i- To understand this 
fact, we first notice that if the bcc operators have the same dimension, hk^i = hp^k,q-i, 
the Liouville momentum is reversed under the replacement {k,l) — )■ {p — k,q — I): 
Pk,i = ~Pp~k,q-i- This sign is normally absorbed by the absolute value in ( |1.2| ). But if 
we take the alternative dressing for B^^q^i, we obtain the relation ^ + f3^q_i = Q- 
This explains the inversion of the Liouville boundary two point function which is 
known to satisfy the reflection relation d{l3\s,t)d{Q — f3\s,t) = 1 pD|. 



In ||3^ was also suggested another symmetry relating the FZZT branes,'^ 

\s; k,l >= —\s — ipq;p — k,l > . (3.19) 



Using the decomposition (|3.11|) , the previous statement is equivalent to 



q-l p-1 

^ \s + imp; 1,1 > = ^ |s + mg; 1, 1 > = . (3.20) 

m=-(g-l),2 n=-(p-l),2 

These linear combinations are exactly those arising in the formal decomposition of 
\s; l,q > and \s;p, 1 >. In the matrix model, we see from ( p.l6|) that the boundary 



one point functions d^Zki vanish whenever k = q or I = p, due to the property 
[q]+ = [p]_ = 0, which supports ( ^.20| ) and the existence of the symmetry ( |3.19| ). 
As a consequence of (|3.19|), the matrix operator Fk{Y) can be equivalently seen as 



creating a boundary with [k, 1) matter BC and Liouville boundary parameter t, or 
(p — k, 1) and t — ipq. 

This dual intepretation of the matrix operator Fk{Y) is also seen at the level of 
one point functions. Let us consider the quotient of the matrix correlators Ok{{y}k) = 
■fjdy (tr log Fk{Y)) and uj{x). Neglecting the non universal terms, we can write 

= (3.21) 



^Note that the Cardy states \k,l > and \p — k,l > have the same boundary entropy up to a sign: 
g(|fc,/>)-(-l)P+9-lg(|p-fc,?>). 
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The RHS contains a product of matter and Liouville contributions. The matter part 
is simply the ratio of matter ipi^i one point functions with {k, 1) and (1, 1) boundary 
conditions, 

(VV)(M) ^ %M)(i.i) ^ . (3.22) 

The second contribution in the RHS of (|3.21|) , the ratio of the resolvents a)(C)/w(0 



is the quotient of the Liouville boundary one point functions with a usual and a 
dual boundary cosmological constants: t = ps/q = IP's. Thus, Ok{{y}k) gives in the 
continuum limit a disc with insertion of a Bi^i boundary operator, with (/c, 1) matter 
BC and t boundary parameter. But the quotient ( p.21j ) allows another interpretation 
when rewritten as 

^44^ = {-lf-'-\p - kU-ir'^-'^ = -i-iy-'-'[p - k]_^, (3.23) 

where we used the property [p — k]^ = (— 1)'^^^[A;]_ and introduced the cosmological 
constants {x, y) with a shifted boundary parameter: x = e^^is—ipq), y = e'^({t—ipq). 
In this setting, Ok{{y}k) is seen as a disc boundary one point function with {p — k, 1) 
matter BC and t — ipq boundary parameter, in agreement with the symmetry (|3.19| ). 
Taking k = 1, we recover the previous interpretation of the resolvant u}{y). 



4. Conclusion 



We briefly summarize the main achievements of this paper and provide some com- 
ments about the further considerations. In the first section, we turned on an new 
light on the construction of the matrix boundaries in the IMM. Starting from a re- 
lation among FZZT branes, we introduced the matrix operator Fi. We showed that 
suitably inserted in a disc correlator, Fi creates a boundary with matter BC (1,/). 
Next, we considered the matrix correlator (9; ^ with boundaries Fi and Fj. and gave 
an explicit proof of the recursion relation relating the critical parts of and 
Ok,i+i ( p.lO| ). This relation provided the identification of Oi^k to a disc with (1,/) 
and (l,fc) matter boundary conditions, and Bi^i^k-i bcc operators inserted. It is 
noted that the construction of all the matrix correlators describing the insertion of 
a different bcc operator between these two matter BC is one of the main remaining 
open questions. We believe that our analysis can easily be extended to more general 
correlators using the linearization provided by the relation ([A.l|) . 

The previous approach was subsequently applied to the 2MM and the boundaries 
with (1, 1) and {p — k, 1) matter boundary conditions were constructed, using respec- 
tively Fi{X) and Fk{Y) matrix operators. As a crosscheck, the matrix correlator 
with Fi{X) and Fk{Y) boundaries was identified to the disc with (1, /) and {p — k, 1) 
matter BC and Bp_k,i bcc operators, in agreement with the CFT predictions. We 



- 13 - 



also clarified the realization of the symmetry {k, I) — ?■ {p — k,q — l) within the matrix 
model, explaining how it is related to the reflection property of the boundary Li- 
ouville two point function. We eventually verified that the symmetry among FZZT 
branes with matter BC {k, I) and {p — k, I) discovered recently in [^| is satisfied by 



the matrix boundaries. This symmetry allows to see the Fk{Y) boundary as wearing 
{k, 1) matter boundary conditions and t Liouville boundary parameter instead of 
{p — k, 1) and t — ipq. The generalization to any {k, I) boundary was also discussed 
but the corresponding explicit matrix correlator needs further research. 

As mentioned in the introduction, to study the boundary RG flows, we have to 
perturb away from the critical points in the space of couplings. These perturbations 
correspond to the introduction of boundary operators that do not change the BC. 
Those remain to be constructed. Moreover, going away from the conformal point is 
a difficult task: The identification of the perturbing coupling of the matrix model 
and the Liouville gravity is far from being obvious. The main complexity lies in the 
appearance of contact terms reflecting the ambiguity of the regularized product of 



operators at the same points These terms lead to a non trivial transformation 
between the bulk matrix and gravity couplings, known as the Belavin-Zamolodchikov 
transform If the transformation has been worked out at the level of the bulk 



one-point function on the disk pH], it still remains to be derived for the boundary 
couplings. 

Let us also mention that the boundaries constructed here for the one and two 
matrix models can also be applied to other models, such as the 0{n) or ADE models. 
It would be very interesting to find out the interpretation of this construction in terms 
of the corresponding statistical models on the fluctuating lattice. The symmetries 
observed above could also be investigated. 
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A. One matrix model 



A.l Technical preliminaries 

In the appendices, we will make an extensive use of the formula 



I 



j = l J j = l J ijLj 

which can be easily proven by noting that the RHS is just the decomposition of the 
LHS as a sum over its poles z = Xj with the proper residue. This formula can be 
applied safely to a matrix as every factor of the product commutes. Note that the 
formula holds for any set of variables Xj,j = 1 ■ ■ - l. 

The equation ([A.l|) can be further exploited by considering its asymptotic when 
z ^ oo. Indeed, expanding the RHS and the LHS separately at infinity, and equating 
both series, one can show that the expressions 

I 

Pa{x) = n ~ ^3^'^ ' a e Z , (A.2) 

are vanishing for < a < / — 1 and are polynomial when a > / — 1. In particular, 
= (-l)'-i. 

A.2 One boundary matrix correlators 

As an example of application for the previous relations, we investigate the matrix 
correlators 

Ol"(W.) = ^(t-^;P^A^'-*) , (A.3) 
The relation ( [A. 1|) allows to expand the ratio l/F^ into 

Of\{^}i)= i: W{x,-x^r^^U.-^M'A. (A.4) 

n=-(/-l),2jyn \ ' 

Then, we recursively make use of the identity 

^ ._M^-^ = M"-^ - M^-^ (A.5) 



Xn- M Xn- M 

in order to decrease the power of the matrix M inside the correlator. The second 
term of the RHS always give a zero contribution because of the vanishing of Pa when 
a <l — 1. We end up with 

l-l 

(^f ^({^}0 = E n (^^- - ^nr'xl-'uj{xn) + (n.u.). (A.6) 

n=-{l-l),2 ji^n 

For the critical choice Xn G {a;};, oj{xn) = (— l)'^^a;(a;) factorizes and the critical 
part of Of"^ is vanishing, unless k = 1: Of\{x}i) = 6ki0o{x) + (n.u.). 
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A. 3 Demonstration of the recursion relation among IMM correlators 

In this appendix, we compute the following two-point function, 




n 



m=-{l-l),2 



M -Xr, 



'■I 
k~l 

n 

n=-(fc-l),2 




We first rewrite Oi u 

clS cl ratio of determinants, 

r -r^ + (n.u.), 



(A.7) 



(Ai 



where A^^^ is given by 
/ 1 



det 



X 



/-I 



1 1 1 

X-(i-\) Vk-i Vk-z 

,2 „ „ 



1 \ 



„A:+«-2 „,fc+Z-2 „,fc+/-2 

i/fc_i yfc-3 



l/-(fe-i) 
V-ik-i) 



y-(fc-i) 



and Ai^fc is the Vandermonde determinant for both variables {x^,?/™} 



(A.9) 



Ai^k = det 



-^1-3 



1 \ 

:fc-i) 

^-{i-i) yk-i Vk-s ' ' ' y-(k-i) 



/ 1 1 ■■■ 1 1 1 

xi-i xi-3 ■ ■ ■ a;_(i_i) yk-i yk-3 ■ ■ ■ y-{k-i) 



^k+l-l „,k+l-l „.k+l-l 

•^-{l-l) i/fc-l yk-3 



(A.IO) 



k+i-i 

y-{k-i)/ 



In writing ( |A.8| ), we assumed that we can replace the resolvents by their critical 
values, up to a non universal polynomial. To show this property, we decompose both 
1/Fi and 1/Fk using the formula ( |A.1|) . The remaining matrix correlator can be 
written using the resolvent W, 

1 1 \ W{yJ-W{xr.) 



Xm - M yn - M 



ym •^n 



(All) 



ym ■^n '^iym •^n) 

In order to replace W by its singular part u, we need to show that the following 
quantities are polynomials, 

vMxu{yh)= E E ^'^'f : r^"^"^ n (^^- - ^-y' n (^^- - ^")" 



m=(-/-l),2 n=-(fc-l),2 



(A12) 
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We notice that the ratio involving the derivative of the potential is actually a poly- 
nomial of Xm and ?/„. It is thus sufficient to show that the quantities 

l-l k-l 
m=(— i-l),2 n=-(fc-l),2 j^n 

are polynomial. In this expression, the two sums can be factorized and we obtain 
the product Pa{{x} i)pb{{y} k) , where the polynomials Pa and pb have been defined in 



We now come back to the expression (|A.8| ) and simplify the determinants using 



the explicit expression of Xm and i/n, Xn = eu cosh (7r6s„) and i/m = eu cosh {nbtr 
with the shifted boundary parameters s„ = s + ibn, = t + ibm. We note that 



x^ 



cosh'=(7r6s„) = {euf2''^^-^^ cosh(fc7r6s„) + 0(x^"^). (A.14) 



In (|A.9|) , we eliminate the terms of 0(x"~^) in ( |A.14| ) by adding other rows in the 
determinant multiplied by some constants. Then, if we subtract the columns of Xm 
and Un from those of Xm+2 and yn+2 respectively and expand the determinant with 
respect to the first and the last rows, we get 

k+l-2 

A,k = coieuY'^^-'^^'^^-'^' [i-ir-My) - i-iy-'^i^)] n 2sin(vr6^j) 

/ sinh(7r6s;_2) ••• sinh(7rfo_(i_2)) sinh(7r64_2) ■■■ sinh(7r6t_(A;_2)) \ 
sinh(27r6si_2) ■■■ sinh(27r6s_(/_2)) sinh(27r6tfc_2) ••• sinh(27r6t_(fc_2)) 

X det 

\sinh(a7r6sz_2) ■■■ sinh(a7r6s_(;_2)) sinh(a7r6tfc-2) ■■■ sinh(a7r6t_(fc_2)) / 

(A.15) 

with Co = {—1)''2~~ ~ and a = + / — 2. At this point, we need the relation 

sin n6' = [/„_i(cos 0) sin ^ where Un{x) is the Chebyshev polynomial of the second 
kind, given by Un{x) = (2x)" + (9(a;"^^). Applying this relation to the determinant 
in ([A.15|), one can again keep only the highest order term in Un{x). Then, we find 



that Ai^k is given by 

A,,fc ={-iy{euf+'-' [{-lf-'uj{y) - {-iy-'u{x)] A,_i (A.16) 

k+l-2 1-2 fc-2 

X JJ 2sin(7r6^j) JJ sinh(7r6s„) JJ sinh(7r6t„), (A. 17) 

j=l m=-(l-2) n=-(k-2) 



Using this formula, one can find that the ratio q''^''"^'^^ is given by a constant indepen- 
dent of s and t, 



Oi,k+2 _ sinh(7r64) sinh(7r&t-fc) A;_i^fc+iA;+2,fc _ [/ + !] + [/]+ 
Oi+2,k sinh(7r6s;) sinh(7rfe_;) Ai+i,fc_iAi,fc+2 [k + l] + [k]+' 
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The relation (|A.18 ) relates Oi^k with k even to Ok+i-2,2 and that with k odd to 
Ok+i-i,i- However, it is easy to show that Om,2 is also proportional to Om+1,1 as 
Cm,2 = [fn]+Om+i,i, so that Oi^k with any k,l is finally reduced to Ok+i-i in the 
continuum limit, 



0,k 



nk+l-2 
j=k 



b1- 



■0 



k+l-1 



(A.19) 



This formula implies the simple recursion relation ( p.lO|) . 



B. Two matrix model 

B.l Derivation of the spectral curve and the mixed trace 

In order to derive the expression of the spectral curve, we use the loop equation 
technique which is based on the invariance of the matrix measure. The first identity 
follows from 



/ 



where we denoted the matrix model action S[X, Y] = tr [V^(X) + U{Y) — XY]. The 
previous equation leads to a relation between the following correlators, 

^Wix) +9 = -j^ (tr ^^U\Y)J , (B.2) 

where we used (tr 1) = A^. 

We need a second relation, obtained from a slightly more complicated expression, 

,XiY tr A ( 1 m^,-.-^,.,) . , (B.3) 

oX \x — X y — Y J 

which in the planar limit, and using (p.2|) , reduces to 
[y - co{x)] I /tr ^ U'{y)-U'{Y) \ ^ u'{y))W{x) - P{x, y) + g , (B.4) 

A' \ X — X y — y I 

introducing the polynomial 

Taking the equation ( [B.4|) at the point y = uj{x), the LHS of (|B.4D disappear and we 
end up with 

E{x, y) = {x- U\y)){y - V\x)) - P{x, y) + g = (B.6) 
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The polynomial E{x, y) is of degree deg V 'm.x and deg U in y, it is called the spectral 
curve. As demonstrated above, it vanishes when y = u){x), or by symmetry x = Cj{y). 



In order to derive the expression of the spectral curve in the double scaling 
limit, we will suppose the simplest choice for the potentials V and U leading to the 
{p, q) critical point: The potential V and U are respectively of degree q and p. As 
consequence, E{x,y) is of degree g in x and p in y, and the equation E{x,y) = 
leads to q solutions Xm at fixed y. Thus, E{x, y) can be written as 



q-l 



E{x, y) = P{y) Ylix- Xm) 



(B.7) 



m=0 



where P{y) is a polynomial in y which will be determined below. Using the parame- 
terization (|3.3|) in the continuum limit, it is easy to see that the condition u{xm) = y 
indeed provide q solutions. 



Xm = 2{euY cosh.{TTSm/q), 



t + 2imp 



(B.8) 



where y{t) = 2{tuY cosh(7rt/]9) and m runs from zero to g — 1. This expression allows 
us to compute explicitely the product 



q-l 



q-l 



\{{x-Xm) = 2'i{tur^\{ 



m=0 



cosh ( — ) — cosh ( — (t + 2imp) 



m=0 '- 

2{euf'' [cosh (its) - cosh (vrt)] 



(B.9) 



where we used the trigonometric formula (2.15) of ||3J]. This expression can be 
written in terms of Chebyshev polynomials of the first kind, Tp(cos^) = cos{p6). We 
thus obtain the following expression for the spectral curve. 



E{x,y)=2{euy'^P{y) 



T 

-J- n 



X 



2(eu)P 



T 



y 



2(eu)i 



(B.IO) 



Since E is of degree p in y, the polynomial P{y) must be a constant. This constant 
can be obtained by comparing the asymptotics of the expression ( [B.IUD and ( |B.(j| ) at 
X — )■ oo. Under a suitable rescaling of the matrices such that V{x) ~ —-x'^ + ■ ■ ■ , we 
can take P{y) = 1, recovering the expression given in 



B.2 Expression for the mixed trace 

To derive the expression of the mixed trace, we start from the following loop equation. 



which leads in the planar limit to 

[x - Q{y)]M{x, y) = -| ^tr ^ ^ ^^'^^l^y^^^ ) " ^(^) • (^.12) 



Plugging the expression (|B.4|) we get 

[x ~ uj{y)){y ~ uj{x)) 
where we recognize the spectral curve in the numerator. 

B.3 Boundary two point function 

In this appendix, we study the continuum limit of the mixed trace matrix correlator 
Mi,k{{x}i, {y}k) defined in ( p.l3| ). Using twice the formula (|A.1| ), and then the 
expression for the mixed trace (|3.8|), one get the double sum 



MMxU{yh) = - J2 1 -f^w^"^ , 

where we recall the notations Xm = x{sm) = 2e^u''^ cosh (nsm/q) and ?/„ = = 
2e'^u'^ cosh {jtn/p) for the cosmological constants with shifted boundary parameters 
Sm = s+imp and t„ = t+inq. The resolvents u and u take the same values for all the 
parameters Xm and ?/„ respectively: co{xm) = {—lY^^u^x) and a)(|/„) = {—l)''^^u{y) 
with X = x{s) and y = y{t). The spectral curve is also independent of m and n, and 
can be factorized out of the sums, 

m=-(/-l),2n=-(fc-l),2 ^ ™ ^ ^ IVyn \ J \ Jj 

(B.15) 

with 

E{x, y) = eV [(-I)p('-^) cosh (tts) - (-l)^^^"!) cosh (vrt)] . (B.16) 
Then we use back the relation (|A.1| ) in order to obtain a factorized form. 



i-i 

MMx}u {y}k) = -E{x, y) n [x„, - {-lf-^uj{y)\ 

m=-(/-l),2 

fe-1 

X n [yn~{-ir'uj{x)]-\ 

n=-(fc-l),2 



(B.17) 
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We now plug in the expression for the cosmological constants and resolvent in terms 
of the boundary parameters, 



Mi,ki{x}i,{y}k) 



i-i 

c. n 

m=-(/-l),2 
fc-1 

X n 

n=-(fc-l),2 



cosh 



TCSr 



- (-D'^-^cosh I — 



fTvt 

V <1 



1, f 

cosh 

V P 



cosh I — 
P 



ns 



-1 



(B.18) 



where we denoted Transforming the difference of hyperbolic 

cosines into a product of hyperbolic sines and reorganizing the indices, one can write 



Mi,ki{x}i,{y}k) 
E{x,y) 



i-i 



Y\ ( sinh 



X 



n=0 



71S± 


+ tlT- 






ns± 




P 


P 



l-l k-lq 

l-lp k-1 
2 q 2~ 



-1 



(B.19) 



with s± = (s ± t)/2 and C2 = 2 ' ^ci. Using the formula (2.15) of in order to 
write the spectral curve as a product of hyperbolic sines. 



E{x, y) = (_i)'?{*^-i)eP%P9 [cosh (vrs - - l)p - iTx{k - l)q) - cosh {-nt)] 

p-i 



± n=0 



h tn- I n 

p p 



l-lp k-1 



2 q 



2 

(B.20) 



we observe that some of the terms simplify. We end up with the following ratio for 
Mi,k{{x}i,{y}k)- 



riLl sinh 



'Ks± q 

h zn- I n 

p p 



l-lp k-1 
2 q 2~ 



^ 11^=0 sinh 



vrs± . p 

h tiT- I m 

q q 



n 



p— fc— 1 

n=0 



sinh 



vrs-t q 

h 111- I n — 

p p \ 2 g 



/ - 1 k-lq 
~^ 2~p 

l-lp p — k — 1 



71S± p 

h ZTi- I m 

q q 



I — 1 p — k — 1 q 
~^ 2 p 



(B.21) 



with the irrelevant constant C4 = 2'^^^^^^^^^\eu)^'^~^^'^~^^ . In the RHS, we introduced 
a shifted boundary parameter, s-t = (s ± t)/2 with t = t — ipq. The expression of 
the Liouville boundary two point function (i(/3|s,t) derived in PD| involve a ratio of 
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double sine functions. But when the Liouville theory is coupled to a minimal model, 
this expression degenerates into a ratio of hyperbolic sine which can be found in 
(appendix A). For a Liouville charge given by ([1.2p, and s and t boundary 
parameters, the formula reads 



Un=l sinh 



d{(3p,a\s,t) = cu 



\Pr,s\y^/2 



n 

^ n^=oSinh 



Trs± . q 

h ivr- I n 

p p 



<T — 1 p p — I 
2 q 2~ 



71S± p 

h ITT- I m 

q q 



a — 1 p — 1 q 
~^ 2~p 



(B.22) 

where c is a constant independent of the cosmological constants. Comparing (|B.21| ) 
and ( |B.22D , one can deduce that in the continuum limit Aik,i is proportional to 
the minimal Liouville gravity two point function of Bp^k,i boundary operators with 
boundary parameters s and t — ipq corresponding respectively to a usual ^(s) and 
dual ~ w) boundary cosmological constants. 

Let us conclude with the following remark: the decomposition of the spectral 
curve can also be made on the dual product of hyperbolic sine. 



E{x,y) = (-l)'?{^-i)+'?-i22('?-i)eP%^'^ JJ JJ sinh 

± m=0 



ZTT- I m — 



l-l k-lq 
~^ 2~p 
(B.23) 



In this case, we obtain another expression for the matrix correlator Aii^k, 



Mi^k{{x}u{y}k) oc W 



n 



r?i=0 



sinh 



vrs-t p 

h m- I m 

q q 



q — I — 1 k — 1 q 
2 2~p 



h tTT- I n — 

p p 



q — I — Ip k — 1 
2 q 2~ 



(B.24) 



which should now be identified with the inverse of ( p.22|) for a B^^q-i boundary 
operator with again s and t — ipq boundary parameters. 
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